Inelastic Scattering in Metal-H2-Metal Junctions 
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We present first-principles calculations of the dl / dV characteristics of an H2 molecule sandwiched 
between Au and Pt electrodes in the presence of electron-phonon interactions. The conductance is 
found to decrease by a few percentage at threshold voltages corresponding to the excitation energy 
of longitudinal vibrations of the H2 molecule. In the case of Pt electrodes, the transverse vibrations 
can mediate transport through otherwise non-transmitting Pt d-channels leading to an increase in 
the differential conductance even though the hydrogen junction is characterized predominately by 
a single almost fully open transport channel. In the case of Au, the transverse modes do not affect 
the dl/dV because the Au d-states are too far below the Fermi level. A simple explanation of the 
first-principles results is given using scattering theory. Finally, we compare and discuss our results 
in relation to experimental data. 

PACS numbers: 73.63.Rt, 72.10.Fk, 85.65.-l-h 
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In recent years it has become possible to measure the 
electrical properties of single molecules captured between 
metallic electrodes [H 0, |j] . Such experiments provide 
a unique opportunity to develop our understanding of 
basic quantum mechanical phenomena at the nanometer 
length scale and at the same time constitute the first 
steps towards molecule-based electronics. [1] 

Interactions between the conduction electrons and the 
molecule's vibrational degrees of freedom is of particu- 
lar interest for the performance of molecular electron- 
ics devices as they determine the local temperature and 
stabilityof the device when subject to an external bias 
voltage[5j. Moreover, inelastic scattering can be used to 
identify the atomic structure of molecular junctions by 
exploiting the sensitiveness of the molecule's vibrational 
frequencies and the electron-phonon interaction to the 
junction geometry. [E S i, S E, Ei El 

Perhaps the simplest molecular junction consists of 
a single hydrogen molecule sandwiched between metal 
electrodes, see Fig. [T] [l^l Shot noise measurements 
on Pt-D2 contacts show that the conductance is car- 
ried predominantly by a single almost fully transparent 
channel[l4], and density functional theory (DFT) calcu- 
lations have shown that this is consistent with a lin- 
ear bridge configuration. [^, 15, [l^, 17 1. An alterna- 



tive configuration where the H2 molecule is dissociated 
in the contact has also been proposed, however, this junc- 
tion yields a conductance larger than IGq [Gq = 2e^ /h 
is the conductance quantum) with contributions from 
three channels jl3|. Inelastic point contact spectroscopy 
provides information about the hydrogen molecule's vi- 
brational frequencies and their variation upon stretch- 
ing. The data obtained from such measurements have 
also been found to be consistent with the linear bridge 
configuration. Q 

The fact that the hydrogen junction supports a sin- 
gle, almost fully open conductance eigenchannel suggests 



that the inelastic scattering processeses should be par- 
ticularly simple to understand. Indeed, consider a junc- 
tion supporting a single scattering channel at the Fermi 
energy with a transmission probability of T = |t(ef^)p. 
At low temperatures the molecule sits in its vibrational 
groundstate and the electron looses the energy Ml to the 
molecule during a scattering event. Assuming a bias volt- 
age eV = /ii — /ifl > Ml an electron incident on the 
molecule from the left with an energy just below /i^, 
must end up in a left moving scattering state after in- 
teracting with the molecule. This follows from energy 
conservation and the Pauli principle. Upon inelastic scat- 
tering, the probability for the electron to enter the right 
electrode is thus changed from T to 7^ = 1 — T. Conse- 
quently, the change in conductance due to the electron- 
phonon interaction should be proportional to 2T — 1, 
i.e. an increase (decrease) in the conductance is expected 
for T < 0.5 (T > 0.5). The same conclusion has been 
reached using more rigorous arguments 1^ 2^, 2l[ and 
has recently been supported by measurements on Pt-H20 
junctions [22l |. 



In this paper we present DFT calculations for the 
dl/dV curves of Pt-H2-Pt and Au-H2-Au junctions in 
the presence of electron-phonon interactions. For both Pt 
and Au electrodes, scattering on the longitudinal modes 
lowers the conductance by a few percentage of Gq in ac- 
cordance with the simple one-channel model discussed 
above. In the case of Pt, the transverse modes can me- 
diate tunneling through the otherwise closed d-channels 
leading to an increase in the conductance of up to 5% 
of Go, demonstrating that the metal-H2-metal junction 
cannot be viewed as a simple one-channel system. For 
Au, the transverse modes have no effect on the conduc- 
tance because only s-states are present at the Fermi level 
and these do not couple via the transverse vibrations. 
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FIG. 1: The supercell use to model the metal-H2-metal junc- 
tion. Only the hydrogen atoms are allowed to vibrate (the 
"dynamic" atoms). This is a good approximation due to the 
large diflterence in mass between Au/Pt and H. The effect of 
the field generated by the vibrating If atoms is taken into ac- 
count inside the indicated inelastic region. The central region, 
C, is coupled to semi-infinite bulk electrodes and periodic 
boundary conditions are imposed in the directions perpendic- 
ular to the contact axis. 



The Hamiltonian of the system is given by 

H — Hc\ + iJph + ^^el-ph, (1) 

where Hc\ is the Hamiltonian of electrons moving in the 
static equilibrimn structure, iJph describes the vibrations 
of the H2 molecule, and i/ci-ph is the interaction between 
the electrons and the vibrating hydrogen atoms. For Hc\ 
we use the Kohn-Sham Hamiltonian. 

Within the harmonic approximation the molecular 
vibrations are described by the Hamiltonian iJph = 
'^xhQ.x{h\h\ + ^) where h\ {b\) creates (destroys) a 
phonon in mode A. The electron-phonon interaction 
takes the form 



H, 



cl-ph 



E E 



(2) 



where the first sum runs over Wannier functions lo- 
cated in the inelastic region, see Fig. [H and the sec- 
ond sum runs over vibrational modes. The electron- 
phonon coupling matrix, M'^, is given by — 
{4>n{r)\W'^ {r)\(l)„i{r)) , where the displacement potential, 
W^{r) — Vus[{R.n}](i")-QA, is the derivative of the effec- 
tive KS potential in the direction defined by eigenmode 
A. In practice is obtained as a finite difference be- 
tween equilibrium Hamiltonians describing the electronic 
system when the hydrogen molecule has been moved in 
the positive and the negative normal direction. 

The current flowing into the molecule (central region 
C) from lead a = L, R is calculated from the formula [28|, 



Tr 



S<(e)G>(£)-I]>(£)G<(£) d£ (3) 



where Gq' is the lesser and greater Green functions of 
the central region evaluated in the presence of coupling 
to leads and the phonons [30| . 

The lesser and greater Green functions are given by 

G^is) = G^(e)[sf (e) + s|(e) + sj,(£)] G'^(e) (4) 
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and 



where G^{e) — [e - 
G^ie) = [G-(£)]t. 

The self-energy originating from the coupling to the 
leads are calculated using standard techniques [27|. For 
the self-energy due to the electron-phonon coupling from 
mode A we use the first Born approximation. 



M^Gf{e±hnx)M^ 



Sph,A(£) = ^[Sp\,(e)-Sp\,(e)] 



e - e' 



(5) 



(6) 



where the last equation follows from the general identity 
Qr _ Qa ^ Q> _ Q< together with the Kramer's Kro- 
nig relation between ImS'' and ReS''. We assume zero 
phonon temperature corresponding to infinite cooling of 
the vibrations, and thus the number of phonons has been 
set to zero. Consequently electrons never interact with 
an excited molecule and therefore can only lose energy to 
the molecule during a scattering event. 

As done often we have omitted the Hartree term in 
the electron-phonon self-energy. 31, sl] The correspond- 
ing energy-independent contribution to the retarded self- 
energy can be understood as a static phonon-induced 
change in the mean-field electronic potential. It is ex- 
pected that this small static potential would be, at least 
partially, screened if included in the DFT self-consistency 
loop. 

The supercell geometry of the considered hydrogen 
contact is shown in Fig. [1] The distance between the 
two electrodes, or equivalently the length of the super- 
cell, has for the case of Pt been chosen to make the 
calculated vibrational frequencies of the H2 molecule 
match the experimental values as close as possible^. 
For the case of Au where less detailed experimental data 
is available, we have chosen the distance by minimiz- 
ing the toal energy. Using the plane-wave pseudopo- 
tential code Dacapo[2^ we have relaxed the surface lay- 
ers, the pyramids and the hydrogen molecule to obtain 
stable junction structures. We used an energy cut-off 
of 25 Ry for the plane-wave explansion, described the 
ion cores by ultrasoft pseudopotentials [24| . and used a 
1x4x4 Monckhorst pack grid for the fc-point sam- 
pling. Exchange and correlation effects were described 
with the PW91 functional. [2^ As a basis for the elec- 
tronic states we use partially occupied maximally local- 
ized Wannier functions [261 which allows for an efficient 
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and accurate calculation of transport properties as de- 
scribed in Ref. 

The vibrational eigenmodes, {Qa}, and corresponding 
frequencies, {^^a}, of the H2 molecule are obtained by 
diagonalizing the dynamical matrix of the system which 
in turn is calculated from the DFT total energies by fi- 
nite differences. Thanks to the large difference in mass 
between the metal and hydrogen atoms, we can calcu- 
late the dynamical matrix for the two H atoms keeping 
all metal atoms fixed. Following this procedure we ob- 
tain a longitudinal stretching mode (Ml), a longitudinal 
center-of-mass mode (M2), as well as two pairwise de- 
generate transverse modes which we refer to as hindered 
rotations (M3) and hindered transverse center of mass 
modes (M4). The modes are sketched in the insets of 
Fig. [2] and the corresponding frequencies are given in the 
caption. 

In Fig. [2] we show the differential conductance calcu- 
lated from Eq. ([3]) including scattering on the different 
vibrational modes separately. To extract the features due 
to the inelastic scattering from those due to elastic scat- 
tering we have subtracted the elastic signal, i.e. we plot 
G{V) = Gfun(^) - Gei(y) + G^i{V = 0), see Ref. M for 
a discussion of this procedure. 

The conductance curves of Fig. [2] present several in- 
teresting features: For both Ft and Au the longitudinal 
modes lead to a decrease in the conductance as expected 
from the one-channel model. It is noticed that the in- 
ternal stretching mode has a much larger impact on the 
electrons than the CM mode. For Au, the transverse 
modes have no effect on the transport, while for Ft they 
lead to an increase in the conductance. Since the junc- 
tion has one fully open channel this seems to conflict with 
the one-channel model which would predict an increase 
only for junctions with conductance < O.SGq. It is noted 
that the differences in the zero-bias conductances are due 
to the tails of the electron-phonon self-energy, which al- 
though centered around the vibrational frequencies also 
have weight at other energies. Before discussing the ori- 
gin of the above mentioned features it is useful to consider 
a simplified description of the scattering process. 

In the following we regard Hei-ph as a perturbation to 
Hq = Hei + i?ph and consider the scattering of a sin- 
gle electron off a molecule in its vibrational groundstate. 
For simplicity we disregard the effect of all the other elec- 
trons (but we do take the Pauli principle into account). 
The scattering states of Hq are conveniently chosen as 
the eigenchannels incident on the molecule from the left, 
ipLpie), or right, ^b.j(£).[33I. The probabihty that an 
electron of energy e injected from the left lead in mode 
p, is transmitted (reflected) upon scattering elastically 
on the central region is denoted by TLp{e) = |tLp(e)p 
{TZLp{e) = |rLp(e)p). Due to the non- mixing property of 
the eigenchannels we have Tap + TZap = 1 for all chan- 
nels p and a = L,R. In terms of the eigenchannels 
the Landauer formula for conductance takes the form 
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FIG. 2: (Color online) Differential conductance of the Pt- 
H2-Pt (full) and Au-H2-Au (dashed) junctions when scat- 
tering on a single vibrational mode is included. The in- 
sets illustrate the vibrational modes together with the sym- 
metry of the corresponding displacement potential VK^(r). 
Frequencies of the H2 vibrational modes (in meV) for Pt: 
hnMi = 190, hnM2 = 171, hnMS = 64, hnM4 = 30 and for 
Au: h^Mi = 249, hnM2 = 141, ftf^MS = 84, h^AM = 37. 



Gci - GoEp'^Lpiep) - GoEg'^Riisp). The state of 
the molecule is specified by the number of phonons in 
each mode, |n). We use the symbol to denote a state 
of the combined electron-molecule system. 

Assume that eV = I-Il ~ fJ'B. > and consider an elec- 
tron incident on the junction from the left in the state 
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V''" = ■0ip(£) with flu < s < fiL and the molecule in 
its vibrational groundstate, |0). According to scatter- 
ing theory, the system ends up in the asymptotic out 
state, ^f™* — S\Lp;0), where S is the scattering opera- 
tor incorporating the effect of -ffci-ph- In the first Born 
approximation we have the transition amplitudes 



{aq;n\S\Lp;0) « (aq; n\Lp; 0) 
'2TriS{Ein - Eout){aq;n\Hci-ph\Lp; 0), 



(7) 



where i?in and -Bout arc the total energies of the combined 
electron phonon system in the in- and out-going states. 
This allows us to express the out state as 



9,A 



where the prime in the sum means that only modes with 
hQx < eV are included. The expansion coefficients are 

c^, = DR,{e - hnx){^Rg{e - hnx)\W^ir)\^,Lpi£)) (9) 

where Dijg(e) is the electronic density of states for 
channel Rq. The normalization constant, Cp = (1 + 
Yl'q xWpql'^)^^^^ 1 has been introduced because the first 
Born approximation is not a unitary approximation to 
S. The fact that only states coming from the right elec- 
trode are included in the sum of Eq. ([8]) is a simple 
consequence of the Pauli principle. 

In the elastic case, an electronic wavepacket con- 
structed from the states ipLp in a narrow interval around 
the energy e, initially located far from the molecule in 
the left lead, will make it to the right lead with probabil- 
ity 7ip(e). In contrast the scattered state (|H]) describes 
a situation where the initial wavepacket makes it to the 
right electrode with probability 



Pp — \Cp\ 



Tlp{£f) 



,A |2 



Ti-Rqiep) 



(10) 



where we have assumed that T and TZ varies little on the 
scale of Kl. The total change in conductance due to the 
inelastic scattering can then be obtained from Landauer's 
formula 



AG = Go 



E 



q.X 



r,^ |2 



'R-Rqisp) - TLpiep) 



ill) 



Apart from the assumptions of instantaneous cooling 
of the phonons and weak electron-phonon interaction, 
which also underlie the first-principles results, Eq. (fTT|) 
was derived in the absence of a Fermi sea. However, as 
we show below, Eq. pT|) provides a simple and physi- 
cally appealing explanation of the first-principles results 

of Fig. da). 

It follows from Eq. pi|l that the change in conductance 
involves all pairs of channels for which the matrix element 
{i^Lp\W^{Y)\'ipiiq) is non-zero for some mode A. Since 



W^{y) extends to the metal atoms binding to H2, any 
scattering state - transmitting or not - with weight on 
these atoms will also contribute in Eq. (fTT|) . 

In the case of Pt, we find at the Fermi level two types 
of eigenchannels with sufficient weight on the hydrogen 
atoms and the contacting Pt atoms that the coupling 
matrix element will be significant. One eigenchannel is 
the almost fully open s-channel and the others have d 
character and very low transmission at £p. Since Au has 
no d-states at the Fermi level, only the s-channel makes 
a contribution in Eq. pT|) . 

For the longitudinal modes, Ml and M2, the symmetry 
of implies that s-s transitions are possible, but not 
s-d transitions {d-d transitions are not excluded by sym- 
metry, but because of the vanishing overlap between ip^d 
and ijjRd)- Since T^.^ — 7^ « — 1 we should expect a drop 
in conductance in agreement with the first-principles cal- 
culations. On the hydrogen molecule, the s-channel has 
mainly character of the II2 anti-bonding orbital. This im- 
plies that the product V'Ls(r)*'!/'/Ss(r) is unchanged upon 
reflection in the plane cutting through the H-H bond per- 
pendicular to the molecular axis. On the other hand the 
potential W^^'^ (r) changes sign upon this reflection. Con- 
sequently, the matrix element {i^ rs\W ps) will be al- 
most zero, and this explains the weak signal observed for 
M2 as compared to Ml. 

The spatial shape of the d-states implies that coupling 
to the s-channel is possible only via the transverse modes 
M3 and M4, see the symmetry of in the insets of 
Fig. [21 Limiting the sums in Eq. pT|) to these two 
relevant states we see that AG becomes proportional to 
TZd — Ts- The increase in conductance found for the trans- 
verse modes in the Pt contact can thus be explained by 
a higher reflection probability of the low-transmitting d- 
channel as compared to the transmission probability of 
the high-transmitting s-channel. We stress that small 
changes in the transmission probabilities for the s- or d- 
channels could change the sign of AG. The symmetry of 
the displacement potential, , for the transverse modes 
prevents coupling between two states with s-symmetry, 
which explains why the transverse modes do not affect 
the conductance of the Au junction. 

We notice that the calculated increase in conductance 
due to the transverse modes is not in agreement with 
the experimental data from inelastic point contact spec- 
troscopy for Pt-Il2-Pt junctions which show a conduc- 
tance decrease. Some of the possible explanations for 
this disagreement is: 

- According to Eq. pT|) . the size (and the sign) of 
AG is determined by the relative magnitude of the 
s- and d-channel transmissions. Even small changes 
here could change the sign of AG. In this sense, 
the fact that we obtain an increase in conductance 
while experimentally a decrease is observed, should 
be viewed as a quantitative rather than a qualita- 
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tive difference. 

- In principle the IBA applies in the limit of 
weak electron-phonon interactions while we obtain 
electron-phonon matrix elements {M in Eq. ([2])) on 
the order of electron volts. On the other hand the 
inelastic features in the dl/dV are a few percentage 
of Go indicating that only a few out of a hundred 
electrons are scattered. Moreover, previous studies 
applying the IBA to gold chains agree nicely with 
experiments^], indicating that the IBA provides 
an accurate description of electron-phonon inter- 
actions in strongly coupled metal-molecule-metal 
junctions. 

- The highly symmetric geometry of the metal-H2- 
metal junction used in this study is an idealised but 
oversimplified model of the real structure. How- 
ever, we have considered other less symmetric con- 
figurations none of which gave rise to a conductance 
decrease for the transverse modes. 

- Inclusion of a finite phonon temperature could af- 
fect the calculated properties. However, as can be 
seen from Eq. (14) of Ref. [35], to lowest order in 
the electron-phonon interaction strength the sign 
of AG cannot change by including heating. 

Despite the differences between the experimental and 
theoretical findings for the phonon induced features in 
the dl/dV, we hesitate to conclude that the linear bridge 
configuration is not the structure observed in the exper- 
iments. The reason is the strong evidence mentioned in 
the introduction which favors the linear bridge combined 
with the small size and high sensitivity of the inelastic 
features. 

In conclusion, we have performed first-principles cal- 
culations for the non-linear dl/dV curves of Pt-H2-Pt 
and Au-H2-Au molecular junctions in the presence of 
electron-phonon interactions. For both metals, the longi- 
tudinal vibrations of the H2 leads to a decrease of the con- 
ductance at bias voltage corresponding to the frequency 
of the vibration, eV — h£l. In the case of Pt electrodes, 
the transverse vibrations induce an increase in conduc- 
tance. This might seem surprising since the hydrogen 
junction supports a single almost fully open transport 
channel and thus, according to the one-channel model, 
inelastic scattering should always lower the conductance. 
On the basis of scattering theory we showed that the in- 
crease is a result of non-transmitting d-channels which 
couple to the transmitting s-channel via the transverse 
modes. This is consistent with the finding that trans- 
verse modes do not affect the conductance in the case of 
Au electrodes. 

We acknowledge support from the Lundbeck Founda- 
tion's Center for Atomic-scale Materials Design and the 
Danish Center for Scientific Computing. 
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